Failure Probabilities and Tough-Brittle Crossover of Heterogeneous Materials with 

Continuous Disorder 
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The failure probabilities or the strength distributions of heterogeneous ID systems with con- 
tinuous local strength distribution and local load sharing have been studied using a simple, exact, 
recursive method. The fracture behavior depends on the local bond-strength distribution, the sys- 
tem size, and the applied stress, and crossovers occur as system size or stress changes. In the 
brittle region, systems with continuous disorders have a failure probability of the modified-Gumbel 
form, similar to that for systems with percolation disorder. The modified-Gumbel form is of special 
significance in weak-stress situations. This new recursive method has also been generalized to cal- 
culate exactly the failure probabilities under various boundary conditions, thereby illustrating the 
important effect of surfaces in the fracture process. 

PACS numbers: 46.30.Nz, 62.20.Mk, 81.40.Np, 83.50.Tq 



I. INTRODUCTION 

Due to their importance in many scientific and engi- 
neering fields, breakdown phenomena in heterogeneous 
materials have been studied extensively in recent yearsE 
It has been well established that the tensile strength of 
heterogeneous materials is dominated by the "weakest" 
regions of the systems, or by the extreme fluctuations 
of the stress fields. Thus the fracture properties of a 
given sample are highly non-intrinsic, non-self-averaging, 
and strongly depend on the specific realization of the 
strengths of individual bonds based on the heterogene- 
ity and microstructure. This dependence of fracture on 
extreme statistics, hence the lack of self-averaging, often 
makes mean- field theories rather dubious. The failure 
probability F n (a) of a heterogeneous system of size n 
under external tensile stress ct or identically the cumu- 
lative distribution function (cdf) of its strength can be 
described by one of the distribution forms of extreme 
value statistics, such as the Weibull distribution 



F n (a) = 1- exp[-n(o-/<T y 



(1) 



or the modified- Gumbel distributiona 

F n (a) = 1 - exphnexpt-A/a 1 ^)] . (2) 

The .modified-Gumbel form was introduced by Duxbury 
et aa in their study of percolation disorder on diluted 
2D fuse networks (square lattice) in which identical lat- 
tice bonds (same conductance and breakdown threshold) 
are present with a probability p. They observed that in 
general the failure of the first few bonds in the diluted 
network worsens the situation fatally and leads to the 
failure of the entire system, or brittle fracture. With 
this result and those of numerical simulations in tough 
and brittle systemscl it was conjectured that the failure 
probability is described for tough fracture by Eq. |l|, and 
for brittle fracture, by Eq. |[ Statistically the degree of 
brittleness of fracture can be characterized by the ratio 



of the size of typical local damage that triggers the sys- 
tem failure to the system size. That is, a fracture is said 
to be brittle if this ratio is small and tough if it is of 
order one. Thus the degree of brittleness of fracture de- 
pends on the disorder type (modeled by thc,distribution 
of the bond strengths) and the system sizeQu. But also 
it depends on the applied stress since the critical dam- 
age shows a dependence on the external stress which we 
shall discuss below. A major difference between perco- 
lation disorder and continuous disorder is that for per- 
colation disorder the failure stress of the system ct;, is 
approximately the fracture stress o\ for the first bond to 
break, while for the continuous disorder one may have 
o~b 3> o~i — > inf{er : G(a) > 0}, where G(ct) is the local 
bond-strength distribution. 

The fracture of materials is typically highly localized in 
the cracks that nucleate and grow due to the local stress 
enhancement at crack tips. In real materials, because of 
the long-range elastic interaction, the stress redistribu- 
tion follows a power law, Act ~ r~ 9 , where Act is the 
stress increase on a bond at distance r from a crack tip. 
The extreme cases of this picture lead to two important 
models which because of their simplicity have received 
considerable attention. First, for g — > one obtains the 
global-load-sharing or equal-load- sharing (ELS) model, by 
which the stress released by cracks is equally shared by all 
the remaining bonds across the sample. This model, ana- 
lytically solvable, is mean field-like, and consequently the 
system-strength has been found to be a Gaussian distri- 
butionB rather than of a type of extreme statistics. Sec- 
ond, for g — > oo one obtains the local-load- sharing (LLS) 
model, by which the stress released by a crack is shared 
only by the intact bonds at the crack tips. This model 
is more realistic in the sense that it does show extreme- 
statistical aspects. In lL\.-an exact, numerical solution 
of this model is availableBQ and it has been used in the 
study of simple 2D systems with linear cracks which can 
be approximately treated as a stack of independent ID 
systems! Here we propose a simple, exact, recursive 
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method which may lead to an ultimate analytical solu- 
tion to the ID problem, and with its help, we obtain 
some insight into the form of the strength distribution in 
more general cases. In addition, by using this powerful 
method, we have re-derived the strength distribution of 
the diluted model analytically. 

Consider an array of n bonds, each of which is as- 
signed a random strength threshold drawn from a con- 
tinuous distribution function G{x). In this article we 
shall use, for example, the Uniform(0,l) distribution: 
G(x) — x, < x < 1 and the Weibull(m) distribu- 
tion: G(x) = 1 — exp(— x m ). We assume the disorder 
is quenched and when the stress on a local bond exceeds 
its assigned strength threshold, it fractures irreversibly 
and quasi-statically. As an external tensile stress a is 
applied, some of the bonds may fracture immediately 
and the stress released will be redistributed among the 
temporarily surviving bonds and this stress enhancement 
might trigger a secondary wave of bond fractures, etc. 
The system fails in one dimension when the last surviv- 
ing bond fails, or in higher dimensions when a spanning 
crack forms. It's easy to see that the failure of a compos- 
ite depends sensitively on how the stress is redistributed 
among the surviving bonds. In the local-load- sharing 
(LLS) model, the stress released by the formation of an 
internal crack is equally shared by the two bonds at the 
tips. However, if the crack is on the end of the sample, 
the simplest method is to assume that there is an intact 
bond outside to bear the stress. This is called an interior 
boundary condition, and is strictly applicable only when 
the sample is embedded in a larger one with intact bonds. 
Thus if a bond has k neighboring failed bonds in the LLS 
model, the stress on this bond is <7fe = (1 + k/2)a, where 
a is the applied stress, and this bond survives with prob- 
ability Wk(c) = 1 — G(crfc). We shall discuss this and 
a number of more complicated and realistic boundary 
conditions, namely, J.he periodic, semi-open, and open 
boundary conditionso. 



II. RECURSIVE SOLUTION OF THE 
ID-LATTICE MODEL WITH LLS 

In this section, we discuss the recursive solution of the 
ID LLS model with interior boundary condition. In order 
to find the failure probability (or strength distribution) 
of the composite given the above, we develop a power- 
ful and simple .recursive method, simpler than the one 
used previously!! To find the failure probability, we need 
to evaluate the sum of the survival probabilities for all 
configurations (except the one with all failed bonds) and 
we define S n ,i to be the sum of the survival probabilities 
for all possible configurations in a sample of size n with 
I fractured bonds on the far right end, so that the failure 
probability can be written F n (a) = 1 — Yli=a Sn,i{o~)- 
To find the recursive relation for S n j, we consider the 
following configuration 
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FIG. 1. S n j(a) as a function of Z, the number of broken 
bonds on the right end. Also shown is Wi(a) vs. I, where 
Wi(a) = 1 - G[(l + l/2)a]. Bond-strength distribution G: 
Weibull(2); system size: n — 100; applied stress a — 0.1. 
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where "1" represents an intact bond and "0" represents 
a fractured bond. Suppose we remove the intact bond 
and the I broken bonds at the rightmost end, we are left 
with a similar configuration, but of smaller size n — l—1, 
therefore the survival probabilities {S n ^i} should satisfy 
the recursive relation 

n-2-2 

S n ,i = F n -i-iW„-iFi + 2^ Sn-l-l,rW r+ lFi , (3) 

r=0 

where as a special case the first term corresponds to the 
surviving configuration having only one intact bond. We 
further define Sk,k to be the failure probability of a such 
system of size k, i.e. Sk,k = F/., then the above results 
can be rewritten in the compact form 

n-l 

S n ,n = 1 — ^ ] S ni l 

1=0 
n-l-1 

S n j = Sn-l^i^Wr+lSlJ , I = 0, • • • , 71 - 1 , 

r=0 

(4) 

with an initial condition Sq.o — 1 which serves as a " seed" 
for the entire calculations. This exact recursion formula 
can be evaluated numerically very efficiently. Despite its 
simplicity, so far we have found that it's hard to analyt- 
ically solve Eq. [I[ Nevertheless, by noticing that both 
Wi and S n ,i decay rapidly with increasing /, and that for 
I small, S n j/ S n -i t i is essentially a constant, we find an 
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excellent approximation to Eq. [|. Fig. 1 shows the I de- 
pendence of S n ,i and Wi for a given system of size n = 100 
with bond strength distribution Weibull(2) under stress 
a = 0.1; it is obvious that these functions decay in some 
exponential manner. Fig. 2 shows the ratio of survival 
probabilities of two systems of sizes n and n — 1 respec- 
tively with I broken bonds on the right end, S n ,i/S n -i,i, 
as a function of I, under the same conditions. The ratio 
is very close to 1 and nearly a constant. 
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G: Weibull(m=2) 
n= 100 
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FIG. 2. 1 — S n ,i{c) I S n -i,i(a) as a function of I. The sur- 
vival probability ratio Sn.i(o') / S n -i.i(a) is very close to 1 
and nearly independent of I. Bond-strength distribution G: 
Weibull(2); system size: n = 100; applied stress a = 0.1. 

From Eq. it is obvious that 



n-l-l 



s, 



n-l-2 



^ S n -l-l,rW r +l — 7T~^ — ^ ^re-l-2,r^r+l i 



r=0 

and also that 



r=0 



-1 n-l-l 



1 — S n>n — ^ ^ S n -l-i :r W r+ lSl : l 



1=0 r=0 

n— 2 g n—l — 2 

i n *n— 1,1 n 
1=0 1 r=0 

,71— 1 • 

So since S nt i/ S n —i t i is nearly independent of /, we obtain 

-2 n-l-2 



i-s n 



Sn,0 
Sn-1,0 



S n .O 



,n— 1 ; 



i=0 r=0 



1 — S„. n — — 1 (1 — S n -i. n -i) + W n -iS n -i tTl -i, 

Jn-1,0 

n=l,2,.-. , (5) 



where the second term W n -iS n -i tn -i is the survival 
probability of the configuration with only one single in- 
tact bond at the left end, and the first and dominant 
term is the total survival probability for all the other 
configurations. 
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FIG. 3. (a) Failure probability as a function of system size 
under different stresses. Open circle: exact numerical data 
by Eq. ^; filled diamond: approximation with Eq. ^. (b) The 
determination of a critical size n c . For n > n c , 1 — S n ,n(&) 
decays exponentially with n, and 1 — S n ,o{&) / S n -i,o((j) be- 
comes essentially a constant. Inset: stress dependence of n c , 

This is an extremely good approximation, which gives 
both the large- and small-n behaviors of S n _ n (a) or F n {a) 
very well, as shown in Fig. 3(a) where the exact results 
of numerically evaluating Eq. ^ are compared with Eq. |^. 
We thus focus on the ratio 5 n ,o/ S n -i,o which gives the 
primary behavior of the fracture for sample of size n. 
First, for sufficiently large n, we notice that in Fig. 3(b), 
S, h o/ Sn-i,o goes to a constant. Thus by neglecting the 
single-bond term in Eq. we obtain that asymptotically 
the survival probability of a system decays exponentially 
with n, which means that a system is always brittle for n 
sufficiently large. Also, Eq. |5| shows its consistency with 
the weak link condition, with which a large system can 
be treated as a collection of smaller sub-systems that are 
weakly related or nearly independent. From Eq. |5|, it can 
be easily shown that 
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(1 >Si,i)(l S n —i in — i) ~ (1 Sjj)(l Sn—j,n—j) ) (6) 

for i, j, n — i, and n — j all much greater than one, 
given S n< o/ S n -i,o goes to a constant. This can be in- 
terpreted as that if a large system is divided into parts, 
then the product of survival probabilities for all parts 
does not depend on how the system is divided, which 
is also guaranteed by the weak-link condition. Second, 
for n small the n dependence of the failure probability 
has a more complicated behavior as shown in Fig. 3(b) 
and the system is tough. This suggests the existence 
of a crossover between different statistics, and allows us 
to specify a tough-to-brittle size scale n c , such that the 
breakdown behavior of a system of size larger than n c 
becomes brittle. In other words, a system of size n 3> n c 
can be treated as a composite of subsystems of size n c , 
which shows brittle-like properties for cracks larger than 
the scale n c . This size scale is analogous to the critical 
size in the nucleation process of a first-order phase tran- 
sition. Given a sufficiently large system, the failure of 
the disordered composite behaves like a first-order phase 
transitionla in that with the nucleation of a critical crack, 
the system fails. It's reasonable to define n c to be the sys- 
tem size beyond which S n $/ S n ~i,o becomes essentially 
a constant [Fig. 3(b)]. The numerical results show that 
n c oc (T _1 independent of the local bond threshold distri- 
bution G(x), but the coefficient does depend on G(x), as 
shown in Fig. 3(b) inset. 

We thus come to a two-scale picture of fracture: lo- 
cally within size n c the fracture is tough and globally the 
fracture is brittle at, scale n c . This confirms the picture 
proposed by Curtintl recently. In his letter, Curtin stud- 
ied planar systems under stresses in the normal direction 
and observed the existence of a critical size n c , according 
to which the system could be characterized as a collection 
of subsystems of size n c . A subsystem in Curtin's picture 
was formulated in an equal load sharing (ELS) model, 
while the whole system became brittle in that failure of 
any single subsystem necessarily led to the breakdown of 
the whole system. The question remaining is to deter- 
mine n c . In our numerically solvable ID model, we find 
a power-law stress dependence of n c oc er" 1 when a is 
small. As we'll describe below, with n c oc er~ Q , the fail- 
ure probability belongs to the " modified-Gumbel" family 
of the extreme statistics. 

For our method we find, and for more general models 
suggest, that for a system of fixed size n with applied 
stress cr, the failure distribution will generally have these 
behavioral regions as follows: 

a) n < n c : This is the tough region, which cor- 
responds to the initial downhill portion of the curves 
in Fig. 3(a). Generally this region is characterized by 
very small stresses (which may be the normal operat- 
ing stresses in practical applications) and smaller system 
sizes. The expansion of Eq. ^ gives the exact failure prob- 
abilities and, in particular it is applicable for very small 
n's. The failure probability is a superposition of a huge 
number of local distributions G(x), and because the fail- 



ure of the composite is path (or ordering) dependent, the 
upper bound of the number can be of order 2 n ~ 1 n\. For- 
tunately, we are only interested in the left tail of the fail- 
ure distribution, which is essential in the extreme (min- 
imum) statistics analysis. Consider a Weibull(m) local 
bond-strength distribution: G(x) = 1 — exp(— x m ), the 
survival probability of a bond with k fractured neighbor 
bonds is W k {a) = exp{-[(l + k/2)a] m }. For the left 
tail of the distribution S n ^ n {o~) where both a and n are 
small, we can make Taylor series expansions of Wk{o~) 
about a = in Eq. ^, and we find the failure probabili- 
ties for samples of different sizes are given by 

51.1 =a m + 0(a 2m ) , 

5 2 . 2 - [-1 + 2 ■ (^)"V m + 0{a 3m ) , 

5*3,3 = [1 - 2 • (^) m - 2 m + 4 • 3 m - (^) 2m ]a 3m 
+ 0{a 4m ) , 

54,4 = [-1 + 2 ■ (^) m + (^) 2m + 2 2m + 2 1+m 
- 2 • 5 m + 2 3 ~ m • 15 m ]cr 4m + 0(a 5m ) , 



For weak stresses, we neglect the higher order terms and 
obtain 

F n (a) = S n , n (a)^c(n,m)a mn 

w 1 -exp(-c(n,m)a mn ) , (7) 

where c(n, m) is a coefficient. It is a Weibull distribu- 
tion with the modulus proportional to the system size 
n. Fig. 4(a) demonstrates the consistency between the 
numerical results and this analytical formula. All lower- 
order terms cancel out in the expansion; physically, this 
can be understood with the fact that the failure of the 
composite requires the fracture of all the bonds, and each 
bond has a contribution of ~ a m to the failure probabil- 
ity. From this point of view, the failure probability for 
the tough region at the low-stress limit is always of this 
form, and the load-sharing rule can only affect the co- 
efficient. The determination of the coefficient c(n,m) is 
hard, but a crude estimate can be made from the n! ways 
of fracture ordering of the n bonds. For small n, due to 
the LLS, we estimate c(n, m) ~ (n!) 7 ™, where < 7 < 1 
is a parameter depending upon m. Thus the failure prob- 
ability for the tough region is of the Weibull form 

F n (a) -l-expHnirV" 1 ") . (8) 

The optimum sample size, corresponding to the minimum 
failure probability as shown in Fig. 1(a), is now easily es- 
timated. By setting F n -\ = F n , we obtain the optimal 
size 

Tlmin ^ 0~ ^ j (9) 
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and the corresponding minimum failure probability 

F mm -expt-ca- 1 ^) . (10) 
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FIG. 4. (a) Failure behavior in the tough region, n < n c . 
Failure probability is a Weibull distribution with modulus 
ran, where n is the system size and ra is the Weibull modulus 
for the local-bond strength distribution, (b) Failure behavior 
in the brittle region, n 2> n c 2> 1. For a very wide range 
of failure probability, from 10 -30 to 0.999, it has a modi- 
fied-Gumbel form, (c) Failure behavior in the super-brittle 
region, n c ~ 1. The strength of the system approaches that 
of the first broken bond. 



For different G(x), we found empirically 7 = 0.59 for Ex- 
ponential, 0.71 for Weibull(2), 0.76 for Weibull(5) and 

0. 89 for Uniform(0,l). 

b) n 3> n c 3> 1 : This is the brittle region, where 
the system is macroscopically brittle but microscopically 
tough. In practice for large system size, most fracture 
events of real materials fall into this region, since the 
other regions are compressed. The discussion on the fail- 
ure probability oLa system with percolation disorder un- 
der a given stressB is still valid in principle for a system 
with continuous disorder when n 3> n c . Roughly speak- 
ing, the failure of the system depends on whether the size 
of the weakest spot exceeds the critical size n c , and the 
probability to find a weak spot of size larger than n c de- 
cays exponentially, so we would still expect a " modified- 
Gumbel" type failure probability for continuous disorder. 
We observe, as shown in Fig. 3(b), for large n and rela- 
tively weak stress, that S n ,o/ S n —\,o varies only with a, 

1. e.. 



Sn-l.O 



(a)^l-exp(clog(a)/a«) 



(11) 



where exponent a = 0.90 ± 0.05 for all local bond distri- 
butions tested, is close to 1, the value expected with the 
load sharing rule LLS. Combining Eq. ||, we get the form 
of the failure probability 



F n (a) = 1 — exp{— anexp[clog(<7)/<7 Q ]} , 



(12) 



where a, c and a are parameters to be determined in data 
fitting. If we use the weak-link hypothesis and treat the 
system as a collection of subsystems of size n c oc cr -1 , 
then the failure probability can be evaluated by using 
F n ~ 1 — (1— F nc ) n /™ c , where F nc is the failure probability 
for the tough region at size n c . We thus obtain a similar 
form to Eq. [l^, which implies that the weak-link hypoth- 
esis is, in principle, correct. With a minor discrepancy in 
a, Eq. |l2| belongs essentially to the "modified-Gumbel" 
family. Thus-for continuous disorder, as well as percola- 
tion disorderu, the failure behavior can still be described 
with the modified-Gumbel form. Fig. 4(b) demonstrates 
that the failure probability fits this form extremely well 
over the very wide range from 10~ 30 to 0.999. 

Although this ID model is only applicable to prob- 
lems of breaking a sheet-shaped object such as a piece of 
paper, we believe and propose that the fracture proper- 
ties of higher-dimensional materials are similar and the 
strength-distribution form of Eq. |l2| is generally applica- 
ble. 

c) n c ~ 1 : This is the super-brittle region of fracture. 
The stress is so strong that critical nuclei exist almost 
everywhere, and thus almost all the bonds fail simulta- 
neously. The failure probability is then simply 



F n (a) = 1 - [1 - G(a)Y 



(13) 



where G(a) is the local bond strength distribution. 
Fig. 4(c) shows that as stress becomes large, the fail- 
ure probability approaches form Eq. 03 asymptotically. 
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It should be noted that the interior boundary condition 
considered here does not guarantee the total stress to be 
conserved, especially for small n, at which a deviation 
from form Eq. ^ can be observed. 

Based on the above discussion, we can roughly sketch 
the phase diagram of the different regions as shown in 
Fig. 5. A tough-brittle-superbrittle crossover occurs as 
external stress increases, while a tough-brittle crossover 
occurs as system size increases. Though the crossover is 
gradual, each region is characterized by a different form 
of failure probability with its own physical interpreta- 
tions. The tough region has a Weibull type failure prob- 
ability, and in this region, the load-sharing rule-seems to 
be not very important. Indeed, in 2D, CurtinB has ap- 
proximated the power-law load-sharing rule with a mean- 
field, global, load-sharing model. The local stress en- 
hancement for this case is quite weak, especially in higher 
dimensions, and the bonds fracture relatively "indepen- 
dently". For the brittle region, where the "modified- 
Gumbel" form applies, the subsystems of size n c are such 
that the fracture of any single subsystem results in the 
global system failure. In the super-brittle region, the sub- 
system size n c further decreases to the order of the lattice 
size due to the external stress, and the strength of the 
system is that of the weakest bond which is distributed 
with the ultimate form of Eq. O. 
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FIG. 5. Different regions of fracture behavior. 
Given the strength distribution (Eq. n2h, the system- 



size dependences of the average breakdown stress (fib) 
and its variance are of great interest. From Eq. [l2[ by 
neglecting the slow-varying factor log(cr) and taking the 
median as the average, we have 



(cr 6 ) ~ (log 71 ) 



-1/c 



(14) 



The variance can be found by considering the change in 
F n {a) by a given factor, by which we obtain 



(15) 



Thus both the strength average and its deviation decay 
logarithmically with increasing system size. 



III. OTHER BOUNDARY CONDITIONS 

In the previous section, we discussed the strength dis- 
tributions of heterogeneous ID fiber-bundles with inte- 
rior boundary conditions, by which each sample (consid- 
ered embedded in a larger one) has both its ends held 
by intact neighboring bonds. The interior boundary con- 
dition gives no surface effects which, practically, are of 
great interest to studies of the strengths of realistic sys- 
tems. Does the presence of a surface affect the strength 
distributions significantly? The answer is yes, and this 
issue has been addressed approximately by Chen and 
Leathu, by using a transfer-matrix method with a "no- 
lone-bond" approximation. In this section, we shall con- 
firm the results by an exact recursive method to calculate 
the failure probabilities under varying boundary condi- 
tions, namely, the periodic, semi-open, and open bound- 
ary conditions. All the boundary conditions concern the 
load-sharing rules at surfaces: interior boundary condi- 
tions (b.c.) promise extra unbreakable bonds outside the 
sample to bear stresses; open b.c. are for isolated sam- 
ples with both ends open; semi-open b.c. have one end 
of the sample open and the other interior; periodic b.c. 
roll up the sample, making a closed circular system with 
no surface at all. 

The configuration of a system sized n must be one of 
the following: 



k I 

(1) 0~2jHJ)~2]0 , k + I = n - 1 (single bond) , 

a 

k I 

(2) 0~^10- ■■010^0 , k + I < n - 2 (double bonds) , 




(3 or more bonds) 



G 



where k = 0, 1, ■ • • , n — 1 and I = 0, 1, • • ■ , n — k — 1 with 
n = 1, 2, • ■ ■. We define W k (a) to be the survival prob- 
ability of a bond with load (1 + k/2)a. Also we define 

\ l (er) to be the survival probability of a system with a 
string of k and I broken bonds on each end respectively, 
and Fn (c) to be the failure probability of the system 
of size n, under certain boundary conditions. The idea 
to obtain the recursion relations is the following: in gen- 
eral, we cut off the • • • 01 and 10 ■ • ■ strings (each has 
only one intact bond) from both ends and are left with a 

Interior: 



Periodic: 



Semi-open: 



Open: 



with the initial condition Fq 1 = 1. Thus the survival and 
failure probabilities of the system are given by 

n — 1 n — k—1 

4° = E E <U . (M) 

fc=0 1=0 

and 

FU = l-S$. (21) 

These recursive relations are inter-connected with each 
other since F«, F^ s \ S*« and appear on the right- 
hand side, so they must be evaluated simultaneously. 
Fig. 6 Shows the exact numerical results of the fracture 
probabilities under various boundary conditions. It is ob- 
served that the following relation holds for most system 
sizes 



truncated sample with interior b.c. whose survival prob- 
ability is available. The survival probabilities of these 
cut-off parts can be obtained by taking into account the 
specific boundary conditions which involve F^ or F^ 
for smaller systems. In doing this, one must be care- 
ful about the load re-assignments on bonds close to the 
surfaces. Special cases, i.e., the single- and double-bond 
systems must be considered separately. Finally the re- 
cursion relations for each kind of boundary conditions 
turn out to be as follows: 



(16) 



(17) 



(18) 



(19) 



> (F„W and F«) > F« , (22) 

while the relation between F^ 1 and F^ depends on sys- 
tem size and applied stress. The differences between the 
fracture probabilities with various boundary conditions 
can be orders of magnitude, especially around the op- 
timal system size. It is thus evident that most of thf 
failures originate from the surfaces, as stated in Ref.Q. 
We also observe a shift of the optimal system size n^ in 
to a smaller value with more boundary influence. How- 
ever, as system size n becomes very large or the stress 
becomes sufficiently strong, these four probability lines 
merge to one, showing an identical asymptotic behavior. 
This means in practice, when n is extremely large, the 
surface effects can be limited, and the interior boundary 
conditions can be used, and the strength distributions 
with all the boundary conditions will have the same gen- 



41/ = +F^W n - l - 2 F^ k _ l _ 2 W n - k - 2 F^ 

n — k — l — 3 n — k — l — r — 3 

+ E E fMS^vM - 

r=0 t=0 

S n,k,l = F k+l W 2k+2l6 k +l,n-l + F^ +l W^F^ fe _ z _ 2 W n-2 
n — k — l — 3 n — k — l—r — 3 

+ E E H%w k+r+l si i i k ^_ 2 ^ t w k+t+l , 

r=0 t=0 

S n,k,l = F^ ) Wk+2lF l S 5 k +l,n-l + F^ ) W n _/_ 2 F J ^ fc _;„ 2 VF„_ fc+ ;_2F/ s) 

n—k—l—3 n — k — l — r — 3 

+ E E ^W^V^^ww/^ , 

r=0 t=0 

S nl,l = Ffc S) W / 2fc+2iF/ S) 4 + ^„_i + F^ ) W / "„ + fc_i_2F^2 fc _;_ 2 W rl _fe + ;_ 2 F ; ( 
ti—k—l—3 n — k — l — r—3 

+ E E 4 s) w2 k+r s i : ) _ k _ l ^ t w t+2l F l {s \ 

r=0 t=0 
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eral form. 
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FIG. 6. Fracture probabilities of systems with different 
boundary conditions under weak applied stress (a = 0.1). 
Inset: Failure probability curves under higher applied stress 
(a = 0.4) showing the merging curves. 



IV. SOME MATHEMATICAL REMARKS ON 
THE STRENGTH DISTRIBUTIONS 

Finally, we discuss the importance of the "modified- 
Gumbel" form of Eqs. || and [l2| We observe that the 
modified-Gumbel form does not itself belong to any of 
the three known distribution forms of extreme statistics, 
but its asymptotic behavior belongs to and is approx- 
imated with one of these three extreme-value-statistic 
distributions. Formally, it's obvious that the modified- 
Gumbel form F(a) = 1 — exp[— an exp(— b/a a )] is the 
cdf of the minimum of an independent samples drawn 
from a parent distribution G(a) = exp(—b/<r a ) (known 
as the Frechet distribution of maxima). So we obtain 
that the asymptotic behavior of F(a) for n very large is 
Gumbel, according to the theory of extreme- value statis- 
tics. However, because of the singularity at a — > 0, the 
most interesting region of fracture, its asymptotic form, 
the Gumbel distribution, which does not have the sin- 
gularity, loses significantly its accuracy. To see this, for 
simplicity, let a, b and a be unity. Then from the the- 
ory of extreme-value statistics, it can be shown that the 
asymptotic form of 

F{a) = 1 - exp[-nexp(-l/cr)] , (23) 

or the domain of attraction as n — > oo is Gumbel type, 
i.e., 

F(a) — » 1 — exp[— exp(log 2 n ■ a - log n)] , (24) 

which serves as approximation of F(a) = 1 — [1 — G(a)) n 
for n sufficiently large. Unfortunately, this ultimate form, 



Eq. [24|, converges very slowly. We here define the conver- 
gence rate of the asymptotic form Eq. |24| as its deviation 
from the exact form F{a) = 1 — [1 — G(a)] n . Thus it can 
be shown that the convergence rate is about ~ 1/ log n 
(as a comparison, the convergence rate can be ~ 1/n 
for a fast-converging system). This makes the ultimate 
Gumbel form of little practical value. One way to deal 
with this difficulty is to make use of a Weibull-type ap- 
proximation with the Weibull modulus dependent on n, 
called the penultimate form, i.e., 



1.0 




0.00 0.01 0.02 0.03 

O 

FIG. 7. The modified-Gumbel distribution function 
F n (a) = 1 — exp[— nexp( — l/c)] and its ultimate asymptote, 
the Gum- 
bel distribution F n (a) = 1 — exp[— exp(log 2 n ■ a — logn)] 
and its penultimate asymptote, the Weibull distribution 
F n (a) = 1 - exp[-(21ogn • a - l) lo s™/ 2 ], with n = 10 20 . 
Inset: A linear-log plot. Observe the differences between the 
distributions for F n (a) small. 

F(a) -> f -exp[-(21ogn-cr- l) lo s"/2] , (25) 

which fits very well in almost the entire range of fail- 
ure probability. The penultimate Weibull form can be 
obtained by comparing three percentiles with those for 
the exact form. In this Weibull form, the modulus as 
a function of system size, is about 26 for n = 10 23 , 
which is quite consistent with experiment results for the 
Weibull exponent in many real samples. In fact, most 
of the extreme- value distributions encountered in practi- 
cal applications can be approximated by the penultimate 
Weibull form very well. This explains why the Weibull 
form has been so widely used in the study of breakdown 
phenomenon. But it's worth pointing out that the penul- 
timate form suggests that a three-parameter fitting in the 
form of 

F{a) = 1 - exp[-(aa - b) m ] (26) 

should be conducted, instead of the two-parameter fitting 
widely used in practice (Eq. |f|). Detailed study shows 
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that despite the fact that the Weibull form is an excel- 
lent approximation for typical average breaking stresses, 
under very small stress there may still be a substantial er- 
ror of up to several orders of magnitude, though the error 
decreases logarithmically with increasing n. As shown 
in Fig 7, the ultimate Gumbel form significantly over- 
estimates the failure probability while the penultimate 
Weibull form underestimates it. In engineering applica- 
tions, tests are usually done under large stresses with 
large failure probabilities, and are extrapolated to the 
weak stresses of normal operating conditions to estimate 
reliability. Thus this error becomes very important, es- 
pecially when n is not very large. The best numerical 
predictor of brittle failure probabilities under very weak 
stresses and hence most operating conditions is thus the 
modified-Gumbel form (Eq. ^ which should be used in 
the reliability analysis of brittle materials. 

V. CONCLUSION 

In summary, we have significantly simplified the re- 
cursive method for exact numerical calculations of the 
strength distribution of the fiber-bundle or ID-lattice 
model with continuous, local, strength distribution and 
local load sharing. We observed qualitatively different 
fracture behavior for systems of different sizes and under 
different stresses. The critical sample size n c is found 
to depend upon the stress applied, a result which should 
generally apply in higher dimensions as well. The re- 
sults here also support the conjecture that tough frac- 
ture can be described with a Weibull form and the brittle 



fracture with a modified-Gumbel form. Tough-to-brittle 
crossovers for fracture are found to occur as system sizes, 
or stresses change. The modified-Gumbel form can be 
approximated by a penultimate Weibull form for inter- 
mediate failure probabilities, but this approximation may 
not be applicable to very small failure probabilities due 
to significant errors. So the modified-Gumbel form is 
of special importance in the reliability analysis of brit- 
tle materials. The important effects of surfaces on the 
strength distributions have also been discussed with a 
similar exact recursion method. 



1 See, for example, Statistical Models for the Fracture of Dis- 
ordered Materials, edited by H.J.Herrmann, and S.Roux (El- 
sevier, Amsterdam, 1990). 

2 P.M.Duxbury, P.L.Leath and Paul D.Beale, Phys. Rev. B 
36, 367 (1987). 

3 B.Kahng, G.G.Batrouni, S.Redner, L.de Arcangelis and 
H.J.Herrmann, Phys. Rev. B 37, 7635 (1988). 

4 Niu-Niu Chen and P.L.Leath, Phys. Rev. B 57, 9319 (1998). 

5 Didier Sornette, J. Phys. A 22, L243 (1989). 

6 D.G.Harlow and S.L.Phoenix, Int. J. Fracture 17, 601 
(1981). 

7 P.M.Duxbury and P.L.Leath, Phys. Rev. Lett. 72, 2805 
(1994). 

8 S.Zapperi, P.Ray, H.E.Stanley and A.Vespignani, Phys. 
Rev. Lett. 78, 1408 (1997). 

9 W.A.Curtin, Phys. Rev. Lett. 80, 1445 (1998). 



9 



